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Abstract 

We consider the stimulated emission of gravitons from an initial state of thermal 
equilibrium, under the action of the cosmic gravitational background field. We find 
that the low-energy graviton spectrum is enhanced if compared with spontaneous 
creation from the vacuum; as a consequence, the scale of inflation must be lowered, 
in order not to exceed the observed CMB quadrupole anisotropy. This effect is 
particularly important for models based on a symmetry-breaking transition which 
require, as initial condition, a state of thermal equilibrium at temperatures of the 
order of the inflation scale. 
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It has recently been argued that the Universe cannot be in a state of "eternal" 
inflation: not only a primordial de Sitter exponential expansion, without a begin- 
ning in time, is impossible [1] , but also there are quantum cosmological arguments 
[2] suggesting for the inflationary phase just the minimum duration required to 
bring in causal contact scales that are not much larger than our present Hubble 
radius Hq^. 

If this is the case, the initial state, whose perturbations are amplified by the 
subsequent inflationary evolution, should be flxed by the dynamics of the pre- 
inflationary era and could differ considerably from the usually assumed vacuum. 
This difference influences the shape of the perturbation spectrum, as first pointed 
out in [3] for the case of tensor perturbations (graviton production), and recently 
discussed in [4] for the scalar perturbation case. 

It is well known, in particular, that any infiationary model based on a tem- 
perature dependent phase transition necessarily requires a homogeneous thermal 
state as initial condition. The existence of initial thermal equilibrium constrains 
the parameters of such models, and a recent discussion [5] suggests that, in their 
context, a sufficient duration of infiation can be arranged only for low enough 
energy scales. The aim of this paper is to point out that also the conventional 
upper bounds on the scale arising from the gravity- wave contribution to the CMB 
quadrupole anisotropy are to be lowered, if the relic graviton background is pro- 
duced from an initial thermal bath rather than from the zero-temperature vacuum. 

Indeed, the infiationary amplification of the vacuum fiuctuations can be prop- 
erly represented, in second quantized formalism, as a process of pair production 
under the action of the external gravitational field [6] , with the produced particles 
necessarily appearing in a final "squeezed vacuum" quantum state [6,7]. If the ini- 
tial vacuum is changed into a state of thermal equilibrium, while the infiationary 
dynamics is left unchanged, the pair production process leads eventually to a ther- 
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mal mixture of'squeezed number" states, instead of the pure squeezed vacuum. 
This modifies the spectral number density of the produced particles, in such a way 
that their energy distribution deviates in general from the flat Harrison-Zeldovich 
spectrum, even in the case of pure de Sitter inflationary dynamics. The isotropy 
properties of the CMB radiation then provide a bound on the inflation scale, which 
depends on the temperature of the initial thermal bath. 

In order to discuss this effect let us recall [6-8] , flrst of all, that the inflation- 
ary particle production can be described in terms of Bogoliubov transformations 
relating, for each mode k, the |in) (fefej^fc) to the |out) (aA;,a^) annihilation and 
creation operators : 

ak = c+(k)bk + c*_{k)bl , al = c_(k)bk + c*_^ik)bl. (1) 

The Bogoliubov coefficients c±{k) depend on the dynamics of the background 
geometry (in particular on the transition from the inflation to the standard decel- 
erated phase), and satisfy | c+ p — | c_ p = 1. By parametrizing c± as: 

c+{k) = coshr/s , c*_ = e^*'^'= sinhr/-, (2) 

the relations (1) can be rewritten as unitary transformations 

ak = T,kbkT,l , al = Sfc^^I (3) 
generated by the squeezing operator *: 

Sfc = exp[(462 - Zfcfef )/2] , Zk = ne^'^" . (4) 

* Note that, although the correct (momentum conserving) transformation 
should involve the two-mode squeezing operator [9], we are considering here the 
simpler one-mode formalism since it gives completely equivalent results for the 
quantities discussed in this paper. 
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The spectral properties of the rehc radiation are usually derived by starting 
from the |in) vacuum state |0), which satisfies bk\0) — 0. The pair production then 
process leads, for each mode, to the squeezed vacuum state \zk) = S/e|0), such 
that ak\zk) = 0. The average particle number can be expressed in terms of the 
squeezing parameter rk as 

Nk= (0|4afc|0) = |c_(/c)|2 = sinhVfc. (5) 

If we start, however, with a number state \nk) in which n particles are already 
present in the given mode, bl.bk\nk) = nk\nk)., we obtain a "squeezed number" [10] 
state \zk,nk) = Efc|nfe), with 

Nk = {nk\alak\nk) = |c_ p(l + nfc) + nfe(l + |c_ p), (6) 

in agreement with the rules of stimulated emission. If we start, more generally, 
with a statistical mixture, 

Pk = J2Pn{k)\nk){nk\ , ^i'n = l, (7) 

n n 

then [11] 

Nk = TtlpkaK) = |c-(k)|2(l + nk) + nk(l + |c_(k)|2), (8) 

where fik = X]n'^-^"(^) initial averaged particle number. The spectral 

number distribution for the (bosonic) massless particles produced from an initial 
thermal bath, at a co- moving temperature T = (3 ^ , is thus obtained by inserting 
into eq.(8) the thermal average number 

fik = {e^'' - (9) 

This result can be checked directly by explicit use of the "squeezed thermal" 
density operator [10], pst = ^Pt^^ , where pt = exp(— (from now on the mode 
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index is to be understood, if not explicitly written; any correlation among modes is 
moreover neglected, as we are treating perturbations in the linear approximation) . 
Indeed, in the convenient representation defined by the "superfluctuant" variable 
X, whose variance Ax is amplified by the squeezing process [9] , pst takes the form: 

Pst{x,x') = {x\pst\x') = 

where cr = exp(— 2r), and n is the thermal average number of eq. (9) (for cr = 1, 
eq. (10) reduces to the usual thermal density matrix in the configuration space 
representation). The computation of Tr(pstb^b) then reproduces exactly eq. (8), 
with |c_p = sinh^ r. 

Note, incidentally, that, unlike N, the entropy growth AS' associated with pair 
production is not affected by finite-temperature corrections to the initial state; for 
a squeezed thermal mixture, AS turns out to be just the same as that obtained 
starting from the vacuum [12,13], namely AS = — Incr = 2r. This can be easily 
verified by computing — Trpgt In for the two-mode generalization of the matrix 
(10), and by subtracting the initial thermal contribution. 

It is also worth noting that, in the real-time formalism of the thermo-field 
dynamics [14], the thermal vacuum is related to the T = vacuum by a kind of 
Bogoliubov transformation, with squeezing parameter = sinh~^[(e^'^ — 1)"-^/^]. 
Such a transformation acts on a "doubled" Hilbert space, obtained by introducing 
fictitious operators associated with each physical operator. In this context, the 
average number of eq. (8) can be recovered, formally, by considering the state 
obtained from the vacuum by the product of two SU{1, 1) Bogoliubov matrices, 
with parameters ri = rt, and r2 = sinh~^ |c_|, provided the relative phase is 
chosen to be i^i — '!?2 = (2m -|- l)7r/4, with m integer. 

We shall now concentrate, in particular, on the stimulated emission of gravi- 
tons from an initial thermal bath, under the action of a changing background 
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geometry which describes the transition from an inflationary phase to the subse- 
quent radiation-dominated and matter-dominated eras. 

The spectral energy density p{u>), which is the variable usually adopted to 
characterize today's distribution of the produced gravitons [8,15,16] is given by 

p{u;)=u;{^)c^u;'N{u;), (11) 

where N is defined in eq. (8), with the corresponding n of eq. (9) (we have 
neglected numerical factors of order unity, and uj is the proper frequency, related 
to the comoving one, k, hy u> = k/a{t), where a is the scale factor of the background 
isotropic metric). 

The Bogoliubov coeflicients c±(a;), connecting the |in) and |out) graviton 
modes for the inflation — > radiation matter transition, have been computed by 
many authors [8,16,17], in the sudden approximation. In such an approximation, 
one ignores the details of the transition among the three different cosmic phases, 
and the particle production is neglected for modes which never "hit" the effective 
potential barrier appearing in the graviton wave equation. As a consequence, the 
Bogoliubov coefficients are not modifled, in this approximation, if an initial phase 
dominated by a thermal radiation bath is inserted before the de Sitter era, since 
the radiation dominated evolution gives no contribution to that potential barrier. 

We then insert the known expression of c_ (a;) in eq. (11) and measure p(a;), as 
usual, in units of critical energy density pc, defining Q{u!) = p{uj)/ pc- By exploiting 
the fact that |c_] > 1 for all the modes undergoing the parametric amplification, 
we finally get (we follow in particular the notations of [17]): 

n{uj, to) ~ GHi^n^{to){ — f~^" coth(^) , a;2 < a; < 
n{u;,to)c^GHi^n^{to){ — f-^''{ — )-^coth{(^) , uo < u < U2. (12) 

UJi U>2 2 

Here /3o"^ = /3"H^o) is the proper temperature of the initial thermal state, adia- 
batically rescaled down to the present observation time to [Po is defined in terms 

5 



of the comoving temperature (3 as /?(to) = /5oi(to)]; ^^7(^0) — 10~^ is the fraction 
of the critical energy density present today in the form of radiation; o; > 1 is a 
coefficient parametrizing (in conformal time) the power-law behaviour of the scale 
factor; Hi = H{ti) is the curvature scale at the time ti marking the end of inflation 
and the beginning of the radiation-dominated era; uq ~ 10~^® Hz is the minimum 
amplified frequency crossing today the Hubble radius Hq^; ioi ~ lO^cJo is the 
frequency corresponding to the matter radiation transition; finally, is the max- 
imum amplified frequency, related to the inflation scale by iO\ ~ V^^^{R\lMpY^'^ 
Hz {Mp is the Planck mass). 

Equation (12) provides the present energy distribution of a gravity- wave back- 
ground of inflationary origin, obtained from a primordial state of thermal equi- 
librium at a proper temperature Z?""*^. In the limit /5 — > 00, and for o; = 1, we 
recover the well-known flat de Sitter spectrum obtained from the vacuum, with the 
usual frequency dependence (~ <^~^) at low energy, due to the radiation- matter 
transition [15-17]. 

The primary effect of the initial finite temperature is to enhance the low- 
frequency graviton production, with respect to the high-frequency sector of the 
spectrum. In this respect, the initial thermal vacuum mimics the effect of putting 
"more power on larger scales" [18], typical of power-law inflation; with the differ- 
ence, however, that the thermal effects are rapidly damped at high a;, for realistic 
values of the rescaled initial temperature /3o. In Fig.l the spectrum has been plot- 
ted for the de Sitter case (a = 1), in such a way as to represent, for various values 
of /?o, the maximum allowed fraction of critical energy density compatible with 
the CMB isotropy. 

The relic graviton spectrum is mainly constrained, indeed, by three kinds of 
direct observations [15]: CMB isotropy, pulsar timing data, and critical density. 
However, as discussed in [15] and [17], the most significant constraint for fiat or 
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decreasing spectra, such as those of eq. (12), turns out to be the isotropy bound im- 
posed at the minimum frequency luq, where it presently imphes ^{loq, to) ~ 10"-'^° 
(we are making here a conservative use of the COBE data [19], as an upper hmit 
on the graviton contributions to the quadrupole anisotropy). Such a condition, 
imposed on eq. (12), provides a bound on the inflation scale Hi, which can be 
conveniently expressed, in terms of the usual spectral index n = 3 — 2a, 



In the limit /9o — > oo this generalizes, to any value of n, the usual isotropy 
constraint on the curvature scale of de Sitter (n = 1) inflation [20], namely 
Hi ~ 10~^Mp (with an improvement of one order of magnitude with respect 
to [15-17,20], due to the use of the more constraining COBE data). The new 
effect, however, is that for flnite initial temperature the maximum allowed scale 
is in general depressed with respect to vacuum production, as illustrated in Fig. 2 
(for (3oUJo <C 1, in particular. Hi scales like f]^^^^'^^^). The inclusion of the ther- 
mal correction is thus expected to modify the existing relations (see e.g. [21]) 
between the power index and the scale of inflation, obained by fltting the observed 
anisotropy on a 10 degree angular scale. 

In particular, for the inflationary models based on a thermal symmetry break- 
ing mechanism, the initial temperature is not independent from the scale 
of inflation itself. Suppose, indeed, that the inflationary phase transition oc- 
curs at an energy scale M, which is the scale at which the time-independent 
vacuum energy becomes dominant (M is related to the curvature scale Hi by 

1/2 

M/Mp = (Hi/Mp) ' , according to the Einstein equations). At earlier times, 
such that P~^{t) > M, the symmetry is restored, and the Universe becomes 
radiation-dominated. The temperature of the initial thermal ensemble, rescaled 
at the beginning (t = U) of inflation, must therefore satisfy /3{ti)M < 1. 

This condition, rescaled down adiabatically at the present time to, provides 




(13) 
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a bound on the spectral parameter /?o, which depends on the duration Z of the 
inflationary phase [Z = a{ti) / a{ti) , where ti marks the end of inflation): 

ff„.„=Mfij'^)m<z(^y . (14) 



^ai ) \M) ~ \Mp ^ 

For Z — > oo this bound is washed out by the inflationary supercooling of the 
original thermal ensemble. However, for models whose parameters are adjusted to 
give just the minimal amount of inflation required to solve the standard problems 

[22], 



Z~Z^,„-e (^^QMGgvJ VlQiOGevJ ^^^^ 
{Trh is the reheating temperature), the bound (14) can be re-expressed in terms 
of the reheating efficiency, Q = T^/^/M, as 



H \ 2 

2 logio(/3o'^o) ~ logio ( ^ j + glogioQ- (16) 



For each given value of n and Q one has then a minimum allowed temperature 
and a maximum allowed scale Hi, which are fixed by the combination of the 
constraints (13) and (16), as illustrated in Fig.3 for three different spectral indices. 
This effect represents a truly "remnant" of the pre-inflationary Universe, in the 
sense of [22]. For Z > Zmim the maximum allowed value of Hi scales up like 
^2/(3+n) g^j^j becomes Z-independent for Puq ~ 1. 

In conclusion, we want to stress that the numerical results presented here 
should be regarded, in many respects, as semi-quantitative results, because of 
the approximations made and of the uncertainty of the experimental data, which 
has not been completely taken into account in our discussion. Nevertheless, we 
believe that already at this qualitative level two important indications emerge 
rather clearly. 

The first is that any fit of the CMB anisotropy in terms of the gravity wave 
background should include a thermal dependence in the spectrum [according to eq. 
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(12)], in order to take into account the possible finite temperature of the initial 
state. 

The second indication is that, even if a thermal phase transition at the GUT 
scale is certainly not ruled out as a possible inflationary mechanism, the effects 
discussed here seem to provide some motivation for investigating the same mecha- 
nism also at lower scales, such as the electroweak one [23] . The constraints on the 
scale imposed by our results seem to challenge, in particular, thermal models of 
inflation obtained in the context of string theory, where the natural scale is very 
near the Planck one. No difficulty seems to arise, on the contrary, in the context 
of a "duality-symmetric" string cosmology [24] where the inflationary phase is the 
dual counterpart of the present decelerating expansion, and starts naturally from 
a flat and cold low-energy vacuum. 
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Figure Captions 



Fig. 1 Maximum allowed spectral energy density [according to eq. (12)] in a relic 
graviton background, produced after a phase of de Sitter inflation, from an 
initial thermal bath at finite temperature T = (3~^ . The rescaled temperature 
I3q^ is measured here in units of ujq — 10~^^ Hz. 

Fig. 2 Maximum allowed inflation scale versus the spectral index n, according to eq. 
(13), for three different values of the initial temperature (3q^ (in units of a;o). 

Fig. 3 Maximum scale and minimum initial temperature for models with the inflation 

rate of eq. (15), and with reheating efficiency Q = 1 and Q = 10""^. The 
allowed region lies below the curves labelled by the spectral index n, and to 
the left of the curves labelled by Q. 
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